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Abstract. Reliable modeling of the atmospheres of cool white dwarfs is crucial for understanding 
the atmospheric evolution of these stars and for accurate white dwarfs cosmochronology. Over the 
last decade ab initio modeling entered many research fields and has been successful in predicting 
properties of various materials under extreme conditions. In many cases the investigated physical 
regimes are difficult or even impossible to access by experimental methods, and first principles 
quantum mechanical calculations are the only tools available for investigation. Using modern 
methods of computational chemistry and physics we investigate the atmospheres of helium-rich, 
old white dwarfs. Such atmospheres reach extreme, fluid like densities (up to grams per cm 3 ) and 
represent an excellent laboratory for high temperature and pressure physics and chemistry. We show 
our results for the stability and opacity of H and C2 in dense helium and the implications of our 
work for understanding cool white dwarfs. 
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INTRODUCTION 

Cool white dwarf stars are among the oldest stars in our Galaxy. Having depleted their 
nuclear energy sources these stellar remnants simply cool-off releasing the internal heat 
to space. The coolest stars have been cooling for billions of years and the properties 
of their atmospheres tells us about the stellar evolution history in our Galaxy. Atmo- 
spheres of cool white dwarfs are not only of interest to astrophysicists but are ex- 
cellent research objects for condensed matter scientists. Density-temperature profiles 
of cool white dwarfs atmospheres are given in figure 1. As helium is less transparent 
than hydrogen, helium-rich atmospheres reach extreme, fluid-like conditions (p of a few 
g/cm 3 ), under which the strong interparticle interactions affect atmospheric chemistry 
and physics. They represent an excellent dense physics laboratory for investigation of 
properties of highly compressed, warm matter in the gas-plasma transition region. 



AB INITIO INVESTIGATIONS 

We aim to improve white dwarf atmosphere models by combining methods of quantum 
mechanics and the theory of fluids [1, 2]. Because of the difficulties in describing the 
properties of fluid-like H/He mixture, until recently the atmosphere models were con- 
structed using the ideal gas approximations. Improvements in computational hardware 
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FIGURE 1. The p - T atmosphere profiles of cool white dwarf atmosphere models of various H/He 
compositions and effective temperatures. The horizontal line indicates the gas-fluid transition region at 
-0.1g/cm 3 . 

and methods have allowed for computation and simulations of matter under extreme con- 
ditions. This allows for derivation of more reliable input physics and chemistry for the 
models. In our investigation we use density functional theory, which is most common 
quantum method used to treat many-particle systems. Most recently, we have investi- 
gated the properties of the negative hydrogen ion and molecular carbon. 

H~ in dense helium 

The stability of the negative hydrogen ion in dense helium is a test case problem 
in modeling of dense, fluid-like H/He mixtures. Bergeron et al. [3] postulated that H 
undergoes pressure ionization in dense helium, which could explain the failure to detect 
the H bound-free opacity in the spectra of cool, helium-rich stars, whose atmospheres 
are expected to be enriched with hydrogen accreted from the ISM. Our calculations show 
a different picture. The negative hydrogen ion becomes more stable in dense helium and 
does not undergo pressure ionization up to a density of 4g/cm 3 . Its ionization energy 
is given in figure 2. It increases with density up to p ~ 2g/cm 3 , and then gradually 
decreases at higher densities. The negative hydrogen ion probably pressure ionizes, 
but at densities much higher than that found in the cool white dwarf atmospheres. The 
increase in the H ionization energy with helium density should produce a shift of the 
H bound-free absorption edge from the near-IR to the shorter wavelengths, making this 
absorption mechanism invisible in the near-IR spectra of these stars. 
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FIGURE 2. The ionization energy of the negative hydrogen ion as a function of helium density. 

C2 in dense helium 

Peculiar DQ stars (DQp) show distorted Swan bands of C2, which we have modeled 
as pressure-shifted by interactions with the surrounding dense helium fluid (Kowalski, 
these proceedings). We then modeled the optical spectrum of DQp star LHS290 [2] 
by shifting the Swan bands by the amount calculated as a function of density in the 
atmosphere (Ar e = 1.6p(eV)). The result is given in figure 3. When we assume a 
pure helium composition the distortion of the Swan spectrum is larger than observed. 
This is because the density in the region of band formation in pure-He/C model is 
~ 0.4 g/cm 3 , which is much larger than the density required to produce the C2 bands 
distortions that would match the optical spectrum of this star (p ~ 0.05 g/cm 3 ). The 
photospheric density is lower when the hydrogen is present in the atmosphere, resulting 
in smaller shifts. A modest amount of hydrogen (H/He=6.75 • 10~ 3 ) is required to 
reproduce the optical spectrum of LHS290. This finding is not unexpected as the derived 
H/He atmospheric composition is consistent with the H/He abundance derived by Kilic 
et al. [4] for the DQp star J1442+4013 of similar effective temperature (r eff = 5737K, 
H/He=2.M0~ 3 ). 
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SUMMARY 

Because of the extreme, fluid-like conditions found inside their atmospheres the cool 
white dwarfs are ideal objects for investigation of properties of warm, dense, H/He- 
rich matter. We apply modern methods of computational quantum mechanics to address 
interesting problems in condensed matter science and astrophysics. We address the prob- 
lem of stability of the negative hydrogen ion and the properties of molecular carbon in 
fluid-like dense helium. At odds with previous speculations, we found that H remains 
stable in dense helium under conditions of cool white dwarfs atmospheres, however its 
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FIGURE 3. Fits to the optical spectrum of the DQp white dwarf LHS290 (thick solid line). All the fits 
are with r e ff = 5800K and log g=8 (cgs). The parameters of the fits are: H/He=0, C/He = 1 .25 • 10~ 7 , and 
Ar e = 1.6p He eV (solid line); H/He=0, C/He = 1.25 • 1(T 7 , AT e = OeV (dotted line); mixed composition 
He/H model, Ar e = 1.6p He eV, C/He = 1(T 6 and H/He = 6.75 • 1(T 3 (dashed line). All synthetic spectra 
well match the entire known spectral energy distribution of the star. The models are described in Kowalski 
[2]. 



bound-free opacity edge shifts to the optical wavelengths, removing this important opac- 
ity source from the near-IR. Investigation of C2 in dense helium shows that in cool white 
dwarfs the optical spectral bands of molecular carbon should show a blueward shift as 
a result of the increase in the electronic transition energy between states involved in the 
C2 Swan transition. This explains the origin of the DQp stars. The observed value of the 
shifts in the DQp stars indicates that the atmospheres of DQp stars are less dense than 
predicted by the pure-He/C models, implying that the atmospheres of these stars are 
probably polluted by hydrogen, which increases the opacity and decreases the density. 
This finding requires confirmation by detailed analysis of the atmospheres of a larger 
sample of DQ/DQp stars. 
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